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Abstract 

Pl^ We present a scheme for laser cooling applicable for an extremely dilute sample of magnetically 

trapped antihydrogen atoms(H). Exploiting and controlling the dynamical coupling between the 
H's motional degrees of freedom in a magnetic trap, three-dimensional cooling can be achieved from 
Doppler cooling on one dimension using the lsi/2 ~ 2^3/2 transition. The lack of three-dimensional 
access to the trapped H and the nearly separable nature of the trapping potential leads to difficulties 
Q_J in cooling. Using realistic models for the spatial variation of the magnetic fields, we find that it 

should be possible to cool the H's to ^ 20 mK even with these constraints. 
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I. INTRODUCTION 

Two years ago the ALPHA collaboration demonstrated trapping of H atoms. []]] This 
result was quickly followed by an improvement in the trapping rate and a measurement of 
the lifetime of H's in the trap; H were shown to be trapped for longer than 15 minutes. {2! 
The long time that the H's remain trapped is important because it opens the possibilities 
to performing measurements that might require several minutes. ALPHA was able to use 
the long trapping time to perform the first measurement of resonant transitions between H 
bound states, jsl Recently, the ATRAP collaboration have also published results where they 
claim to have trapped H.j4| 

One of the difficulties in measuring the spectroscopic transitions in H is that the trapped 
H's are currently at relatively high energies which will lead to line broadening from the 
Doppler effect and from Zeeman shifts. To date, the trapped H are formed using three body 
recombination js, 6| when an antiproton (p) is inside of a positron (e"*") plasma. Because three 



body recombination is a relatively slow process at the densities and temperatures in the H 
traps compared to coUisional slowingj^, Isj, it is estimated that the p will approximately be in 
thermal equilibrium with the e"*" plasma before recombination occurs. Typical temperatures 
reported for the e"*" plasma have been a couple lO's of Kelvin. Since the magnetic trap depth 
for H is only ~ 1/2 K, the H's have an energy distribution that extends from to the trap 
depth. This was measured in Ref. ^ where the distribution of H annihilations as a function 
of time or as a function of axial position of the trap matched that for the expected H energy 
distribution. 

Currently, trapped H atoms have energies up to ~ 500 mK, and occupy a large volume of 
order (2 cm)^ x 30 cm. They are in a strongly non- uniform magnetic field, varying by order 
1 t|i|, from the center to the walls of the trap. While a number of important experiments 
have been performed or planned in such a trap, laser cooling of H, if achieved, will provide 
a major experimental advantage. 

Laser cooling will create a cold, and spatially localized sample of antimatter atoms. Local- 
ized atoms will be much less susceptible to the Zeeman effect, currently a dominant limitation 
for microwave spectroscopy [sl, and one of the limitations for future laser spectroscopy jol. 
The lower velocities of cooled atoms will reduce 2nd order Doppler broadening for ls-2s 
two photon spectroscopy. Importantly, laser cooling will greatly increase the sensitivity for 
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observing gravitational interaction of antimatter (see e.g. Ref. 

Only one experiment so far has reported laser cooling of atomic hydrogen, nearly 20 years 
ago In that experiment, ~ 10^^ atomic hydrogen at 80 mK, pre-cooled via evaporative 
cooling, was laser cooled to 8 mK in 15 min. While it is theoretically possible to cool the 
H's in a similar manner, H laser cooling presents considerable experimental challenge for 
several reasons: (1) generation of coherent radiation at 121.6 nm remains technologically 
difficult, due to the lack of convenient lasers and nonlinear crystals at these wavelengths. (2) 
experimental requirements for H trapping allow only limited optical access to the trapped 
atoms. (3) because of the very low densities of H, three-dimensional cooling assisted by 
collisionally mixing the degrees of freedom (needed for laser cooling of atomic hydrogen llj ) 
is prohibitive. (4) because of the large Zeeman effects, only a small portion of the trapped 
atoms resonantly interact with photons. 



Several proposals exist for overcoming some of these challenges 12|-|l4l| , but none has 
been experimentally realized. Reference [l^ demonstrated laser cooling of magnetically 
trapped Na atoms to ~ 2 mK using a one-dimensional optical molasses. This experiment 
showed that it is possible to obtain substantial cooling even with the severe restriction to 
one laser. In this paper, we will investigate whether this simpler scheme for laser cooling 
of H will work, and show, via detailed numerical calculations for ALPHA-type apparatus, 
that three-dimension cooling to ~ 20 mK should be possible within realistic experimental 
and technological constraints. 

In our scheme, the Doppler cooling[l6| will drive the lsi/2 — 2p3/2 transition with the 
light being linearly polarized perpendicular to the quantization axis in order to drive the 
m = l/2tom = 3/2 transition. By driving this transition, the light scattering does not 
lead to a spin- flip which would cause the atom to be ejected from the trap[ll|. 

While a powerful narrow-line cw Lyman-alpha laser could eventually offer advantages in 
laser cooling, development of such sources remain considerable challenge jl?!. fisl ]. In this 
work, we consider the use of a modern pulsed Lyman-alpha source 191] , whose time-averaged 
(as well as instantaneous) power is much greater than cw. Because Lyman-alpha generation 
requires highly non-linear processes, a pulsed scheme offers overall better cooling efficiencies, 
as long as the transition per pulse is not saturated. 

A key feature of the present scheme is the exploitation, and the control of the dynami- 
cal coupling between the z- and xy-degrees of freedom. This will allow three dimensional 
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cooling in an ALPHA-type apparatus in which optical access is currently limited to one di- 
mension due to constraints for efficient magnetic trapping as well as high sensitivity particle 
detection. [10 1 



While the trapping fields are clearly non-separable in x, y, z near the walls of the trap, 
this is not true near the trap center. The effective potential energy for the H's near the trap 



center approximately has the form Vi{z) + V2(r) where r = x^ ^ and z is along the trap 
axis. Since the light is nearly parallel to the trap axis, this nearly separable potential allows 
fast cooling of the 2;-motion but can lead to heating in the xy-coordinates. 

However, there is some small coupling between the z-motion and the xy-motions. This 
coupling is the conduit through which we can achieve cooling in all directions. We enhance 
this coupling by the use of non-harmonic magnetic fields in both xy- and z- directions, in 
contrast to standard harmonic magnetic traps. In the xy-directions, the effective potential 
is given by ~ r^, while in the 2;-direction we use total of five solenoidal coils to produce the 
nonlinearity. Three dimensional cooling is possible when the time between photon scatterings 
is comparable to or longer than the mixing time between all of the degrees of freedom. This 
leads to non-trivial behavior of the final temperature on the laser power. 

In this paper, we will present results on many of the important aspects for laser cooling 
in this constrained geometry. We have investigated the time dependence of the cooling, the 
energy distribution versus detuning, the optimum detuning, etc. In the results section, we 
give physical reasons for the difficulties in trying to laser cool H. 



II. NUMERICAL METHOD 

The basic physical situation is that the H's classically move through the trap. When the 
laser is on and the atom is within the waist of the field, the atom can scatter a photon. The 
atom receives two momentum kicks for each time a photon is scattered: when the photon is 
absorbed the atom is kicked in the z-direction and when the photon is emitted the atom is 
kicked in a random direction. The size of each momentum kick is /i/A ^ 5.45 x 10^^''kg m/s. 
This corresponds to a velocity kick to an H of ~ 3.3 m/s. The H's cool when the total 
momentum kick is opposite the momentum of the H. In this section, we describe the 
computational techniques we used to model this process. 
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A. Classical motion 



The H's move through a magnetic trap where their de Broghe wavelength is much smaller 
than the trap dimensions. This means we can solve for their motion using classical forces. 
The potential energy for the center of mass motion is equal to U = —pi ■ B where /i is the 
magnetic moment of the H in the Is state; is approximately the magnetic moment of the 
e^. Since the precession frequency is much higher than other motional frequency scales, the 
angle between pi and B is an adiabatic invariant. This means the orientation of the positron 
spin with respect to the magnetic field does not change. Thus, the trapped H's experience 
a potential U = fiB where the fi is the magnitude of the magnetic moment and B = \B\. 

To compute the force, we need to obtain 

F = -VU = -iJ,VB (1) 

where we need to compute the gradient of the magnitude of the magnetic field. Because the 
magnetic field is a very complicated function of the coordinates, we computed it numerically 
using a central two point difference: 

F^ = -[U{x + dx/2,y,z) - U{x - dx/2,y, z)]/dx (2) 

and similar operations for Fy and F^. We used dx = dy = dz = R x 10~^ where R is 
the radius of the trap ~ 2.2 cm. Although this seems a crude approximation, the error is 
actually quite small. The error term in the gradient is {dx'^ /A8)U"'] since U'" ~ U' /R"^, this 
approximation gives a relative error of ~ 10^^^ which is comparable to the round-off error 
in this approximation. 

One of the big problems in the calculation is that we need to solve for the H motion for 
lOO's of seconds. We need to be careful that there is no energy drift in the calculation which 
would either give an unphysical cooling (which would lead to overly optimistic results) 
or an unphysical heating (which would lead to a suppression of the laser cooling). We 
found that the adaptive step-size Runge-Kutta alg orithm that worked well for the shorter 



times needed to model the results in Refs. \1 



201 1 was not accurate enough for the present 



calculations unless we used very small time steps. We found that the fourth order symplectic 
integrator 



21 



231 worked well for this calculation. As is usual with symplectic integrators, 



we found that the energy varied during the calculation but the variation remained within a 
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small energy region; thus, there was no energy drift at long times. We used a time step of 
20 fis in our calculations. 

As with our previous investigations, we launched the H's within an ellipsoidal region with 
a fiat spatial distribution. The ellipsoid had a scale length of 0.8 mm in the xy-coordinates 
and scale length of 8 mm in the 2;-coordinate. The initial velocity distribution was chosen 
from a thermal distribution with a temperature ~ 50 K except where explicitly stated 
otherwise. Since the trap depth is only ~ 0.5 K, our effective velocity distribution is flat 
in velocity space within a sphere delimited by the trap depth. Before turning on the laser 
pulses, we had the H's move through the trap for 2 s plus a random time between and 0.2 s 
to model the fact that there is a delay between the H formation and manipulations done to 
them. This time delay allows the H to reach somewhat random regions of phase space. 

The trapping field is generated using an octupole field to provide radial confinement and 
mirror coils to provide axial confinement. We used the approximations in Appendix A of 
Ref. {20! for the fields from mirror coils and the octupole field. Instead of two mirror coils, 
we used 5 coils in order to mimic the more complicated magnetic geometry in the ALPHA- 
II trap. All coils are approximated as two loops with the approximate vector potential of 
Eq. (A.2) of Ref. |2oj. All coils have a radius a = 45.238 mm and A = 0.90230. The two 
end mirror coils have a loop separation of 8.425 mm and a center-to-center separation of 
274 mm. The other 3 coils are equally spaced between the two end mirror coils. All 3 of 
these coils have a loop separation of 8.083 mm. 

The purpose of the extra 3 coils is to provide a flatter magnetic field in the center of the 
trap. Without the extra coils, the B-field near the center of the trap is quadratic in z — Zmid- 
One coil at the center can cancel the quadratic dependence and give a B-field that has a 
quartic dependence. The coils at the 1/4 and 3/4 position together with the middle coil can 
cancel both the quadratic and quartic dependence. This will give a B-field proportional to 
{z — ZmidT '^^ar the center of the trap. 

We use 3 different currents through these 5 coils: the end mirror coils have the same 
current and the coils at the 1/4 and 3/4 position have the same current. In most of the 
calculations, we chose currents to give the flattest possible B-field. For our coil parameters, 
this is achieved with 606 A in the end coils, -57.8 A in the coils at the 1/4 and 3/4 position, 
and -2.5 A in the middle coil. We did perform calculations with somewhat different currents; 
these calculations will be discussed in the Results section. 
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Figure 1 shows a slice through the magnetic field given in terms of the Zeeman shift 
in mK units. The nearly rectangular nature of the contours highlight the nearly separable 
nature of the potential energy experienced by the H's at low energy. 

B. Light scattering 

To perform realistic calculations, we need to use parameters for the 121.6 nm light that are 

n 

within current technical capabilities. We used parameters suggested to us by T. Momose.|19| 
We assumed a laser with a 10 Hz repetition rate. We assumed the laser pulse would be on 
for a short time, ~ 10 ns; this time is so short that we consider both the absorption and re- 
emission to happen instantaneously. We used a laser linewidth (FWHM) of 100 MHz. The 
total energy in one laser pulse was taken to be 0.1 /iJ. Calculations were mostly performed 
with the laser propagating exactly along the 2;-axis and the laser was assumed to be linearly 
polarized with a waist radius of 10 mm. Some calculations were performed at different 
directions for laser propagation and will be discussed in the Results section. 

We used a semiclassical treatment of laser cooling to model the interaction between H's 
and the laser field. If the atom is within the laser waist, the probability for absorbing a 
photon during one of the laser pulses is 



where Tgp is the spontaneous decay rate of the 2p state (27r99.7 MHz = 626 MHz), Tias is 
the laser line width (27r 100 MHz, FWHM), Au is the detuning of the laser (combination of 
laser detuning, Doppler shift, and Zeeman shift), Eias is the energy in the laser pulse, and 
w is the radius of the laser waist. We will measure the laser detuning from the transition at 
the minimum B-field, Bmin, and denote it by Auq. The Doppler shift is —uJis,2p'Vz/c where Vz 
is the z-component of the H velocity and u}is^2p = {E2p — Eis)/h is the transition frequency 
between the Is and 2p states. The Zeeman shift is —^b{B — Bmin)/^ where is the Bohr 
magneton because the m = 3/2 upper state shifts more strongly in the magnetic field than 
the m = 1/2 ground state. Thus, the detuning in Eq. ([3]) is given by 



P 




-p V -^V ' ids J j V / -^Las 



(Tsp + Tias)/{2n) Bias 



(3) 



Au = Auq - Uis,2pVz/c- 11b{B - B^in)/h 



(4) 



where B is evaluated at the position of the H. 
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The algorithm to incorporate the photon scattering worked in the following way. We 
stepped the H using the symplectic time step of 20 fis. A laser pulse is sent through the 
trap every 5000**^ step. Nothing happens if the H is outside the waist. If the H is within 
the waist, we use a random number generator and compare to the probability to scatter a 
photon in Eq. If the random number is smaller than this value, then the atom's velocity 
is modified by the two kicks. The two kicks give a change in velocity: v ^ v + zAv + OAv 
where Av (~ 3.3 m/s) is the photon momentum divided by the H mass and the random 
direction is opposite the direction of the photon emission. The vector z> is randomly chosen 
from the photon emission distribution for a circularly polarized state. 

The cooling discussed below should be compared to the cooling in a three-dimensional (B- 
field free) optical molasses. The lowest average atom energy is when the laser detuning is set 
to Auo = —(Tsp + Tias)/2 and gives an average energy of Eav = 3h(Tsp+Tias)/4:. For our laser 
parameters, the Eav/ks — 7.2 mK. The recoil energy is E^ec/kB = MAv^j {2kB) — 0.64 mK. 

III. RESULTS 

There are two important issues that need to be addressed: what is the best laser detuning 
and how long is needed to get substantial cooling. The optimum detuning depends on the 
laser power for this system because the laser only directly cools one direction which is 
nearly separable from the other two directions. For our laser parameters, we found that the 
detuning should be substantially shifted from the B-field free optical molasses. We found 
optimal cooling with the field free detuning when we performed calculations where the laser 
could cross the trap at a large angle so that there was substantial components in z- and x- 
or y-directions, but this will not be an option for the ALPHA experiment. 

We can obtain an overview of the cooling through the time dependence of the temperature 
of the trapped H's. Figure 2 shows the average energy of the laser cooled H's as a function 
of the time that the laser is on. The atoms start with the distribution that is trapped from 
a 54 K H distribution. The figure shows the results for several possible detunings of the 
laser. The solid, dotted, dashed, and dot-dash lines are for the trapping potential in Fig. 1. 
Discussed below, the dash-dot-dot-dot line is the best cooling that could be obtained when 
the trapping potential is only from the two mirror coils at ±137 mm. The value of the 
detunings are given in terms of the optimal value of —{Tias + Tsp)/2 when B = 0. 
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This figure shows that the most rapid coohng occurs at early times and the coohng 
rate has substantially slowed near the final times. A pessimistic interpretation is that this 
suggests it will take quite a long time to reach the asymptotic temperature. As will be seen 
in the figures below, a large fraction of the H's seem to cool to energies under 100 mK within 
200 s while some H's seem to remain at a few 100 mK. The hotter atoms have substantially 
fewer photon scatterings compared to the colder atoms. Thus, the colder atoms will interact 
more with a laser in a spectroscopy experiment which will lead to a somewhat cooler effective 
temperature. 

We also performed calculations for magnetic fields similar to the geometry of the original 
ALPHA experiments. An example is the dash-dot-dot-dot line in Fig. 2. In this case, the 
three middle mirror coils are off and the trapping potential has a quadratic dependence on 
the 2;-coordinate. We found the cooling to be much worse in this case for two reasons. The 
quadratic dependence of the magnetic field on z means the Zeeman shift is substantial for 
a larger region of space; this leads to less photon scattering and, hence, a smaller cooling 
rate. The other reason is that the quadratic magnetic field gave even less coupling between 
the z- and a;|/-directions; this led to a longer mixing time and, hence, a smaller cooling rate. 
As with the flat 5-ficld case treated in this paper, we found much improved cooling if the 
laser is not constrained to be nearly along the 2;-axis. 

Figure 3 shows the time dependence of the energy distribution for different time windows 
over a 200 second duration. The data in Fig. 3 is for when the laser detuning is 4x the 
field free optimal value of — (Fjas + ^sv)l% since we are using Vias — F^p, this means the 
laser detuning is ~ — 4F5p. An important point is that the distribution is clearly giving 
more H's at lower energy as the atoms are cooled longer. This means that substantial laser 
cooling is possible for this trap geometry over a time scale where H's can be trapped. The 
average energy of the H's has decreased from ~ 330 mK to ~ 110 mK over this time. A 
small fraction, ~ 2 — 3%, of the atoms are lost in the coohng process. These are all atoms 
that were barely trapped and the first couple of photon scatterings gave heating instead of 
cooling due to the random nature of the laser cooling. 

The efficiency of the cooling is naturally of interest. The average number of photon 
scatterings by an H is ~ 5 for each 40 s time window. This gives ~ 25 scatterings for the 
full 200 s of our simulation. There are 2000 laser pulses during this period which means 
that there is shghtly better than a 1% chance for scattering a photon in a laser pulse. For 
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TABLE I: Average number of scattered photons, Np, during 200 s cooling and average final energy 
Ef as a function of laser detuning, Aloq in units of —{Fias + Tsp)/'2. For this case, the starting 
distribution was that expected from the ALPHA experiment. 







Ej (uiK) 


2x 


~ 50 


200 


3x 


~ 35 


135 


4x 


~ 25 


115 


5x 


~ 15 


140 



comparison purposes, the speed of an H with 330 mK of kinetic energy is approximately 
75 m/s and the velocity kick from a photon absorption or emission is approximately 3.3 m/s. 
Thus, each photon scattered provides a substantial amount of cooling for this detuning and 
for this duration of cooling. 

Figure 4 shows the energy distribution of the H's during the final 40 s time bin as a func- 
tion of the detuning of the laser in units of the optimal B-field free detuning. Plots are shown 
for detunings of 2x, 3x, 4x, and 5x. All cases started with the same energy distribution of 
H's. It is clear that the detuning leads to strongly differing energy distributions. Prom this 
figure, it appears that the best detuning is 3x the field free value because the peak of the 
distribution is at the lowest energy. However, the average final energy is actually lowest for 
the 4x detuning as seen in the Table I. The average initial energy is approximately 340 mK 
for our simulated trap. 

Prom Table I, we can gain some insight into how the coohng process works for this H 
trap. Smaller detuning leads to more photon scattering but the photon is more likely to be 
scattered by atoms with smaller \vz\ because a smaller Doppler shift brings the photon into 
resonance. The average change in energy during a single scattering is AE = M(f^Af + Af ^). 
The l^zl can be small when the H is cold or when the atom is moving nearly perpendicular 
to the ^-direction. Thus, the small detuning leads to a lot of scattering without much 
energy removed during the scattering event. Since there is relatively little time between 
each scattering event, the H does not have sufficient time to mix the motion in x, y, z. But 
having too large a detuning leads to a different cooling problem. Simply put, there are too 
few photons scattered to give effective cooling during the 200 s simulation time. 



10 



Although Fig. 4 seems to clearly favor the 3x detuning, Table I gives similar average final 
energies for 3x, 4x, and 5x detuning. This is because the larger detuning more strongly 
cools the higher energy part of the distribution. This suggests that the optimal strategy 
might be to change the frequency of the laser so that we start with large detuning at early 
times and change the frequency to smaller detuning at late times. We found that this did 
provide more cooling over the fixed frequency calculation but it was not a qualitative change. 
In our calculations, we changed the frequency linearly with time. For 200 s, the best case 
we tested started with 6x detuning initially and finished with 3x detuning. However, the 
average final energy was 90 mK. Thus, changing the detuning is probably worth doing only 
if it is experimentally easy. 

One question we wanted to address is what is the final energy distribution if one could 
cool for very long times. To address this, we started with cold thermal distributions and 
allowed them to interact with the laser pulses for 200 s. We found that atoms with an initial 
temperature of 50 mK cooled to ~ 30 mK in 200 s with either the 3x or the 4x detuning. 
Therefore, we are presenting results when the atoms start with a thermal distribution at 
temperature 30 mK. 

Figure 5 is the same as Fig. 3 except for the initial velocity distribution. In Fig. 5, we 
start with a thermal distribution with a temperature of 30 mK {E/kB = 45 mK) in order 
to probe what are the lowest temperatures that are achievable. For this case, the energy 
distribution has settled into its final value at late times. We find that the energy distribution 
is well approximated by a thermal distribution at a temperature of 20 mK {E/ks — 30 mK) 
which is shown in the inset. We note that the usual optical molasses temperature would give 
a temperature of 4.8 mK for optimal detuning. Thus, the final temperature for this trap 
and detuning is only a factor of ~ 4 higher than could be achieved with a 3-dimensional 
molasses. We do not have calculations for the time required to reach the final distribution 
when starting from the high temperature case of Fig. 3. However, our estimates indicate 
that it should be less than the 1000 s trapping seen in Ref. ^. 

Figure 6 is the same as Fig. 4 except for the initial velocity distribution. The 3x and 
4x detuning clearly give better final distributions. Only the 2x detuning is still evolving 
at the final time. That case had the distribution evolving to higher energy at late times. 
As in Fig. 4, the 3x detuning has the peak at slightly lower energy than the 4x detuning. 
However, the average final energy is essentially the same for the two detunings. 
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TABLE II: Average number of scattered photons, Np, during 200 s cooling and average final energy 
Ef as a function of laser detuning, Aloq in units of —{Fias + Tsp)/'2. For this case, the starting 
distribution was thermal at 30 mK. 





Np 


Ef (mK) 


2x 


~ 85 


46 


3x 


~ 45 


33 


4x 


~ 25 


32 


5x 


~ 15 


36 



Table II shows similar data to that in Table I except the starting energy is a thermal 
distribution with a temperature of 30 mK. Except for the 2x detuning case, the average 
energy was approximately constant after 120 s which means Table II gives the asymptotic 
average energy for all detunings except the 2x case. Since almost all of the H's are initially 
trapped, the average initial energy is 45 mK. The trends present in Table I are reflected in the 
data of Table II as well. We think these results are quite encouraging for laser cooling since 
the average final energy for a three dimensional optical molasses is ~ 15 mK for the laser 
parameters we used in this simulation. The average energy for the 2 x detuning case was still 
increasing at the final time; we also found that the average energy for the 2x detuning was 
increasing at the final time when we started with a 50 mK thermal distribution. The final 
average energy in that case was ~ 60 mK which means the 2x detuning has an asymptotic 
average energy which is more than two times higher that for 3x or 4x detuning case. 

We have performed calculations for other magnetic field geometries although we do not 
present their details. We tried to increase the cooling rate by increasing the coupling of the 
motion in the z-direction with the xy-directions. We increased the coupling by deliberately 
making a small, non-flat potential in the central region. We observed an increase in the 
cooling rate when wc increased the current in the central coil to make a potential hill at the 
center of the trap. We observed a decrease in the cooling rate when we made a potential 
dip at the center of the trap by decreasing the current in the central coil. However, the 
increase/decrease of the cooling rate was only apparent when starting with low energy H's 
because the size of the perturbations we tested was at the ~ 10 mK scale. 
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IV. CONCLUSIONS 



We have performed calculations related to prospects for laser cooling trapped H atoms. 
Although the standard methods of laser coohng will work equally well for H, the experimental 
restrictions related to access to the atoms, the large magnetic fields present in the traps, and 
the small wavelength of the light require accurate modeling to address how much cooling is 
possible in practice. Our calculations use accurate magnetic field geometries and realistic 
laser parameters. We have found that an asymptotic temperature only a factor of 2 higher 
than for a three-dimensional optical molasses is possible. We only simulated the case where 
the laser light was impinging on the H's along the trap axis; much better cooling is possible 
if the laser direction is substantially away from 0° or 90° relative to the trap axis. Small 
angles did not have a large effect. 

In this paper, we only presented results for nearly fiat magnetic fields because this is 
clearly the geometry that will be desired for spectroscopic measurements. If it is possible 
to use strongly different S-fields and then have them morph to the flat geometry, then 
much lower temperatures could be possible. For example, one might set up a much tighter 
flat region in z using the 5 mirror coils. After coohng, the S-field could be changed to 
that in Fig. 1 which would give adiabatic cooling due to the expansion of the trap region. 
Since these possibilities seem likely to be highly machine dependent, we will save these more 
complicated situations for when the experiments are attempted. 
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ported by the US National Science Foundation and the Natural Sciences and Engineering 
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by the Auburn University Office of the Vice President for Research through the Undergrad- 
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FIG. 1: The potential energy experienced by an H due to the spatially varying magnetic field. 
Every contour represents a change in energy of 10 mK xkn- 
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FIG. 2: Evolution of the average energy of the H's as a function of the detuning when the magnetic 
potential is that shown in Fig. 1. The linetypes are for when the laser is detuned 2x the B-field free 
optimum value (solid line), 3x detuned (dotted line), 4x detuned (dashed line), and 5x (dash-dot) 
line. The dash-dot-dot-dot line is for the optimum detuning when only the two end magnets are 
energized which leads to a trapping potential quadratic in z — Zmid- 
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FIG. 3: Energy distribution of the H's when the laser is detuned 4x the B-field free optimum 
value and the trapping potential is that in Fig. 1. The different curves correspond to different time 
windows: 0-40 s is sohd, 40-80 s is dotted, 80-120 s is dashed, 120-160 s is dash-dot, and 160-200 s 
is dash-dot-dot-dot. 




FIG. 4: Energy distribution of the H's for the time window of 160-200 s when the laser is detuned 
2x the B-field free optimum value (solid line), 3x detuned (dotted line), 4x detuned (dashed line), 
and 5x (dash-dot) line. 
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FIG. 5: Same as Fig. 3 except that the initial energy distribution was a thermal distribution at 

30 mK (current actual plots arc for 50 mK until the data comes in). The inset shows the energy 
distribution in the final window (solid line) and a 20 mK thermal distribution (dotted line). 
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FIG. 6: Same as Fig. 4 except that the initial energy distribution was a thermal distribution at 
30 mK. 
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